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We report updated branching fraction measurements of tire color-suppressed decays B'^ — > D^tt", 
£)*%", D°7?, D'^r), D^uj, D'^w, D°r]', and_D*°?7'. We measure the branching fractions (xlO^*): 
B{B° -^ DV) = 2.69_± 0.09 ± 0.13, B{B° -^ D*°7r°) = 3.05 ± 0.14 ± 0.28, B{B° -> D^i]) = 
2.53 ± 0.09 ± 0.11, B{B° -^ D'^rf) = 2J9 ± 0.14 ± 0.23, B{B° -> D°uj) = 2.57_± 0.11 ± 0.14, 
B{B° -^ D*°u)) = 4.55 ± 0.24 ± 0.39, B{B° -^ D'^ri') = 1.48 ± 0.13 ± 0.07, and S(B° -^ D*°rj') = 
1.49 ± 0.22 ± 0.15. We also present the first measurement of the longitudinal polarization fraction 
of the decay channel D'^ui, /i = (66.5 ± 4.7 ± 1.5)%. In the above, the first uncertainty is statistical 
and the second is systematic. The results are based on a sample of (454 ± 5) x 10 BB pairs 
collected at the T(4S) resonance, with the BABAR detector at the PEP-II storage rings at SLAC. 
The measurements are the most precise determinations of these quantities from a single experiment. 
They are compared to theoretical predictions obtained by factorization. Soft CoUinear Effective 
Theory (SCET) and perturbative QCD (pQCD). We find that the presence of final state interactions 
is favored and the measurements are in better agreement with SCET than with pQCD. 

PACS numbers: 13.25.Hw, 12.15.Hh, ll.30.Er 



I. INTRODUCTION 



Weak decays of hadrons provide direct access to the 
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lia^Y interaction scattering in the final state [1[ (Final State 
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UK must be well understood. The two-body hadronic B de- 
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«,, .,, TT ■ -i- J- n ■ n ■ Ti 1 is a light meson, are of great help in studying strong- 

'Also with Universita di bassari, bassari, Italy . ° . , . ,,,..,„ ,„ , 



interaction physics related to the confinement of quarks 
and gluons in hadrons. 

The decays B — >■ D^*^h can proceed through the emis- 
sion of a W^ boson following three possible diagrams: 



external, internal (see Fig. [T|), or by a W^ boson ex- 
change whose contribution to the decay rate is expected 
to be much smaller than the external and internal am- 
plitudes 0. The neutral B° — >• D'-*^^hP decays proceed 
through the internal diagrams |3| . Since mesons are color 
singlet objects, the quarks from the W^ decay are con- 
strained to have the anti-color of the spectator quark, 
which induces a suppression of internal diagrams. For 
this reason, internal diagrams are called color- suppressed 
and external ones are called color- allowed. 



(a) 



B^ 




gO w 
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c p(*)0 (b) 
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FIG. 1: External (a) and internal (b) tree diagrams for B" 
D^*'h decays. 



We already discussed factorization models ^ ^ 

previous publication [3] . Within that approach the non- 
factorizable interactions in the final state by soft glu- 
ons are neglected. The matrix element in the effective 
weak Hamiltonian of the decay B — >■ D^*^h is then fac- 
torized into a product of asymptotic states. Factorization 
appears to be successful in the description of the color- 
allowed decays Q- 

The color-suppressed & — >■ c decays i?" — > D'*^°7r" were 
first observed by Belle H] and CLEO ^ with 23.1 x 10^ 
and 9.67 x 10^ BB pairs respectively. Belle has also ob- 
served the decays D^rj and D^uj and put upper limits on 
the branching fraction (B) of D*°-q and D*°uj [i]._Thc 
branching fraction of the color-suppressed decays S" — >■ 
^(*)o^o n(*)o^^ £,(*)Oa;,_and D"rf were measured by 
BaBAR [3| with 88x10^ BB pairs and an upper limit was 
set on B{B° -^ D*"ri'). Belle updated with 152x10^ BB 
pairs the measurement of B{B^ — > Z?(*)°/i*'), h° = 7r°, 77, 
w [ni, and 77' 01 and studied the decays B° -)> D°p° 
with 388 X 10^ BB pairs [IJ]. In an alternative approach, 
BaBAR IJI used the charmless neutral B to K^tt^tt^ 
Dalitz plot analysis with 232x10® BB pairs, and found 
S(i?° — > ZJ^tt") to be in excellent agreement with earlier 
experimental results. BaBAR has also performed a pre- 
liminary Dalitz-plot anal ysis of the mode B^ — > ZJ^tt+tt^ 
with 471 xlO® BB pairs |il. 

Many of these branching fraction measurements are 
significantly larger than predictions obtained within the 
factorization approximation [3|, |l6[. But, while the ini- 
tial various experimental results demonstrated overall 



good consistency, the most recent measurements pub- 
lished by Belle [ll|, [l2l have moved on average towards 
lower B for the color-suppressed B^ — > D'-*'^h'^ de- 
cays, closer to factorization predictions. However, it 
has been demonstrated [l7| that non-factorizable con- 
tributions are mostly dominant for the color-suppressed 
charmed S" -^ D^tt'^ decay and therefore cannot be ne- 
glected. 

Stronger experimental constraints are therefore needed 
to distinguish between the different models of the color- 
suppressed dynamics like pQCD [perturhative QCD) Il8l . 
|l9| or SCET {Soft Collinear Effective Theory) |20l422| . 
Finally, we emphasize the need for accurate measure- 
ments of hadronic color-suppressed B° — >• D^*^^h^ decays 
to constrain the theoretical predictions on B^,js decays 
to D^*^P and D^*^P states, where P is a light pseu- 
doscalar meson such as a pion or a kaon [2J] . Using flavor 
S't/(3) symmetry, the comparison of Bd and Bs decays 
offers new possibilities to determine the decay constant 
ratio fs/fd ,2s|- These decays are and will be employed 
to extract the CKM-angle 7 and other angles [2J|, es- 
pecially in the context of the _B-physics program at the 
LHC. 

This paper reports improved branching fraction mea- 
surements of eight color-suppressed decays i?° -> 
£,(*)o^o^ £)(*)%, D(*)Ow, and D(*^°tj' with 454x10® BB 
pairs and presents for the first time the measurement 
of the longitudinal polarization for the decay mode to 
two vector mesons B'^ — )> D*'^ui, which also constrains 
QCD models and challenges Heavy Quark Effective The- 
ory {UqET) (see Sec. EO. 



II. THE BABAR DETECTOR AND DATA 
SAMPLE 



The data used in this analysis were collected with the 
BaBAR detector at the PEP-II asymmetric-energy e+e" 
storage rings operating at SLAC. The BaBAR detector is 
described in detail in Ref. [231 . Charged particle tracks 
arc reconstructed using a five-layer silicon vertex tracker 
(SVT) and a 40-laycr drift chamber (DCH) immersed in 
a 1.5 T magnetic field. Tracks are identified as pions 
or kaons (particle identification or FID) based on like- 
lihoods constructed from energy loss measurements in 
the SVT and the DCH and from Cherenkov radiation 
angles measured in the detector of internally reflected 
Cherenkov light (DIRC) . Photons are reconstructed from 
showers measured in the CsI(Tl) crystal electromagnetic 
calorimeter (EMC). Muon and neutral hadron identifi- 
cation is performed with the instrumented fiux return 
(IFR). 

The results presented are based on a data sample of 
an integrated luminosity of 413 fb~ recorded from 1999 
to 2007 at the T(45) resonance with a e~'"e~ center-of- 
mass (CM) energy of 10.58 GeV, corresponding to (454± 
5) X 10® BB pairs. The equal production rate of B^B^ 
and B~^ B~ at that resonance is assumed in this paper. 



as suggested by the Particle Data Group (PDG) y^. A 
data sample o£41.2 fb"^ with a CM energy of 10.54 GeV, 
below the BB threshold, is used to study background 
contributions from continuum events e~^e~ -^ qq {q = u, 
d, s, c). We call that latter dataset off-peak events in 
what follows. 

Samples of simulated Monte Carlo (MC) events are 
used to determine signal and background characteristics, 
to optimize selection criteria and to evaluate efficien- 
cies. Simulated events e+e" -^ T(4S') -^ B+B' , B°B°, 
e"'"e~ —^qqiq = u, d, s) and e+e^ — > cc are generated 
with EvtGen [231, which interfaces to Pythia [2^ and 
Jetset [111. Separate samples of exclusive i?° — )• D^*^°h^ 
decays are generated to study the signal features and 
to quantify the signal selection efficiencies. We also 
use high statistics control samples of exclusive decays 
B^ -^ £)(*)°7r^ and D''*^"p^ for specific selection and 
background studies. We study these control samples both 
in data and in the MC, using the same selection criteria. 
All MC samples include simulation of the BaBAR detector 
response generated through Geant4 [S^. The equivalent 
integrated luminosity of the MC samples is about three 
times that of the data for BB, one time for e+e" — > qq 
{q = u, d, s) and twice for e+e~ — >■ cc respectively. The 
equivalent integrated luminosities of the exclusive B de- 
cay mode simulations range from 50 to 2500 times the 
dataset. 



III. ANALYSIS METHOD 

A. General considerations 

The color-suppressed B'^ meson decay modes are re- 
constructed from Z?(*)o meson candidates that are com- 
bined with light neutral-meson candidates hP (tt'', 77, w, 
and 77'). The £)(*)" and hP mesons are detected in various 
possible decay channels. In total, we consider 72 different 
^0 _^ £)(*)0/j0 j^gg^y modes. 

We perform a blind analysis: the optimization of the 
various event selections, the background characteriza- 
tions and rejections, the efficiency calculations, and most 
of the systematic uncertainty computations are based on 
studies done with MC simulations, data sidebands, or 
data control samples. The fits to data, including the 
various signal regions, are only performed after all analy- 
sis procedures are fixed and systematic uncertainties are 
studied. 

Intermediate particles of the decays B^ — > D^*'''^hP 
are reconstructed by combining tracks and/or photons 
for the decay channels with the highest decay rate and 
detection efficiency. Vertex constraints are applied to 
charged daughter particles before computing their invari- 
ant masses. At each step in the decay chain we require 
that the candidate mesons have masses consistent with 
their assumed particle type. If daughter particles are 
produced in the decay of a parent meson with a nat- 
ural width that is small relative to the reconstructed 



width, we constrain the mass of this meson to its nominal 
value, except for the lu and the p° [2g. The 5° mass is 
computed using the constraint of the beam energy (see 
Sec. IIIICTj) . This fitting technique improves the resolu- 
tion of the energy and the momentum of the B'^ candi- 
dates as they are calculated from improved energies and 
momenta of the D(*)° and hP . 

Charged particle tracks are reconstructed from mea- 
surements in the SVT and/or the DCH, and they are 
assigned various particle identification probabilities by 
the FID algorithms. Extrapolated tracks must be in the 
vicinity of the e'^e~ interaction point, i.e. within 1.5 cm 
in the plane transverse to the beam axis and 2.5 cm along 
the beam axis. The charged tracks used for the recon- 
struction of 77 — >■ 7r"'"7r~7r*^ and 77' — > ii^ tt^ r}{-^ 77) must 
in addition have a transverse momentum pT larger than 
100 MeV/c and at least 12 hits in the DCH. When a 
PID positive identification is required for a track, the 
track polar angle 6 must be in the DIRC fiducial region 
25.78° < 6 < 146.10°. Photons are defined as single 
clusters of energy deposition in the EMC crystals not 
matched to a track, and with shower lateral shape con- 
sistent with photons. Because of the high machine back- 
ground in the very forward part of the EMC, we reject 
photons detected in the region 9 < 21.19°. We assume 
that the production point of the photons is the recon- 
structed primary vertex of each e^e~ collision. 

The selections applied to each meson (tt^, ry, uj, rj', I?°, 
and D*") are optimized by maximizing the figure of merit 
S/y/S + B, where S is the number of signal and B is the 
number of background events. The numbers S and B 
are computed from simulations, and the branching ratios 
used to evaluate S are the present world average values 
of color-suppressed decay modes [2al . Each particle mass 
distribution is fitted with a set of Gaussian functions or 
a so-called modified Novosibirsk empirical function |31| . 
which is composed of a Gaussian-like peaking part with 
two tails at low and high values. Particle candidates are 
then required to have a mass within ±2.5 a around the 
fitted mass central value, where a is the resolution of 
the mass distribution obtained by the fit. For the decays 
D^ — > K^Ti^Ti^ and D*^ — > D^^/, the lower bound is 
extended to 3 cr because of the photon energy losses in 
front of and between the EMC crystals, which makes the 
mass distribution asymmetric with a tail at low values. 



B. Selection of intermediate particles 

1. tt" selection 

The t:^ mesons are reconstructed from photon pairs. 
Each photon energy E^ must be greater than 85 MeV 
for TT*^ produced directly from i?" decays, and greater 
than 60 MeV for 7r° from 77, uj, or D'^ meson decays. 
Slow neutral pions originating from D*^ —J- D^ir^ decays 
must satisfy E^/ > 30 MeV. The 7r° reconstructed mass 
resolution ranges 6.5 — 7.0 MeV/c^ for tt*^ from rj, lu, and 



D^ mesons decays, and 7.0 — 7.5 MeV/c'^ for vr" produced 
in Z?*° or B" decays. 



2. Tj selection 

The 77 mesons are reconstructed in the 77 and tt+tt^tt*' 
decay modes, accounting for about 62% of the total decay 
rate [2^], and may originate from S" -^ L*'*'"?] or if ^• 
7r+7r~77 decays. 

The 7? — > 77 candidates are reconstructed by com- 
bining two photons that satisfy £'7 > 200 MeV for B^ 
daughters and E^ > 180 MeV for 77' daughters. As pho- 
tons originating from high momentum tt" mesons may 
fake a 77 ^ 77 signal, a veto is applied. The 77 — >■ 77 
candidate is rejected if either photon combined with any 
other photon in the event with E^ > 200 MeV has an in- 
variant mass between 115 and 150 MeV/c^. Such a veto 
retains 93% of the signal while reducing the background 
of fake 77 mesons candidates by a factor of two. The reso- 
lution of the 77 — > 77 mass distribution is approximately 
15 MeV/c^, dominated by the resolution on the photon 
energy measurement in the EMC. 

For 77 candidates reconstructed in the decay channel 
TT^TT^TT^, the TT^ IS required to satisfy the conditions de- 
scribed in Sec. HUB II The mass resolution is about 
3 MeV/c^, which is better than for the mode 77 —J- 77, 
thanks to the relatively better resolution of the track- 
ing system and the various vertex and mass constraints 
applied to the rj and tt*^ candidates. 



3. uj selection 

The w mesons are reconstructed in the 7r+7r~7r'' de- 
cay mode. This mode accounts for approximately 89% 
of the total decay rate. The tt" is required to satisfy 
the conditions described in Sec. HUB II and the trans- 
verse momenta of the charged pions must be greater than 
200 MeV/c. The natural width of the uj mass distribu- 
tion r = 8.49 MeV [2y] is comparable to the experimen- 
tal resolution cr ~ 7 MeV/c^, therefore the w mass is not 
constrained to its nominal value. We define a total width 
(Ttot = ■\/cr^ + r^/c* ~ 11 MeV/c^ and require the uj can- 
didates to satisfy \muj — Hm^ | < 2.5 atot (where /i,„^ is 
the mean of the w mass distribution) . 



4-. p^ selection 

The p^ mesons originate from 7/' -^ p'^j and are recon- 
structed in the tt+tt" decay mode. The charged tracks 
must satisfy pt(7''^) > 100 MeV/c, where pt is the trans- 
verse component of the momentum with respect to the 
beam axis. We define the helicity angle 9pO as the an- 
gle between the direction of the momentum of one of the 
two pions and that of the 77' both evaluated in the p" 
center-of-mass frame. Because the p*^ is a vector meson. 



the angular distribution is proportional to sin OpO for sig- 
nal, and is flat for background. The p° candidates with 
I cos0po| > 0.73 are rejected. Due to the large p° natural 
width r = 149.1 MeV ^, the mass of the p° candidate 
must lie within 160 MeV/c^ around the nominal mass 
value and no mass constraint is applied. 



5. 77 selection 

The 77' mesons are reconstructed in the tt^tt^ij{^ 77) 
and p°7 decay modes. These modes account for approx- 
imately 46.3% of the total decay rate. 

Only the 77 — ^ 77 sub-mode is used in the 7r"'"7r~77 re- 
construction due to its higher efficiency. The selection 
is described in Sec. HUB 21 For candidates reconstructed 
in the p'^j decay channel we select p° candidates as de- 
scribed in Sec. IIII B 41 and the photons must have an en- 
ergy larger than 200 MeV. As photons coming from tt'^ 
decays may fake signal, a veto as described in Sec. HUB 21 
is applied. The 77' mass resolution is about 3 MeV/c^ for 
7r+7r~77 and 8 MeV/c^ for p°7. 



6. K% selection 



The A'g mesons are reconstructed through their decay 
to two charged pions (7r~7r"'') which must originate from 
a common vertex, with a x^ probability of the vertex 
fit that must be larger than 0.1%. We define the flight 
significance as the ratio L/ctl, where L is the K^ flight 
length in the plane transverse to the beam axis and a^ 
is the resolution on L determined from the vertex fit. 
The combinatorial background is rejected by requiring a 
flight significance larger than 5. The reconstructed K^ 
mass resolution is about 2 MeV/c^ for a core Gaussian 
part corresponding to about 70% of the candidates and 
5 MeV/c^ for the remaining part, depending on the trans- 
verse position of the K^ decay within the tracking system 
(SVT or DCH). 



7. -D" selection 

The D" mesons are reconstructed in the i^~7r+, 
K~'k~^it'^, K~Tr'^Tr~'iT'^ , and /v^tt+tt^ decay modes. 
These modes account for about 29% of the total de- 
cay rate. All D'^ candidates must satisfy p*{D^) > 
1.1 GeV/c, where p* refers to the value of the momentum 
computed in the T(45') rest frame. That requirement is 
loose enough so that various sources of background can 
populate the sidebands of the signal region. 

For the decay modes reconstructed only with tracks, 
we require that the charged pions originating from the 
D° candidates fulfiU ptItt^) > 400 McV/c for K-Tr+, 
PriT^"^) > 100 MeV/c for K~TT+n~TT+, and priTr^) > 
120 MeV/cfor K^7r-T:+. 



The charged tracks must originate from a common ver- 
tex, therefore the x^ probabihty of the vertex fit must be 
larger than 0.1% for the decay channel K^n^ and larger 
than 0.5% for the other modes with more abundant back- 
ground. Because of the increasing level of background 
present for the various decay modes, the kaon candidates 
must satisfy from looser to tighter PID criteria for the 
modes K~Tr^, iir~7r+7r~7r+, and K~'!T~^Tr^ respectively. 
For 7^g7r+7r~, the K'^ candidates must satisfy the selec- 
tion criteria described in Sec. IIIIB6I 

For the decay D'^ — > A'^tt+tt" the combinatorial 
background can significantly be reduced by using the 
parametrization of the K~Tr~^Tr'^ Dalitz-plot distribution 
as provided by the Fermilab E691 experiment [32|. This 
distribution is dominated by the two K* resonances 
{K*° -> K-TT+ and K*' -^ K'-k^) and by the p+{-k+'k°) 
resonance. Therefore we select only D^ candidates that 
fall in the enhanced region of the Dalitz plot as deter- 
mined by the above parametrization. The tt*' must satisfy 
the selections described in Sec. IIIIB II 

The reconstructed D^ mass resolution is about 5, 5.5, 
6.5, and 11 MeV/c^ for the decay modes K^tt^tt^tt^ , 
iiTgTr+Tr^, Ar^7r+, and K'n^Ti^ respectively. 



D*^ selection 



the B candidate where the B energy is set to the beam 
energy in the CM frame: 



m-ES 



En 



- \pi 



(1) 



The variable AE is the energy difference between the 
reconstructed B energy and the beam energy in the CM 
frame: 



AE = El,„ + El - yi/2. 



(2) 



where y^ is the e+e~ center-of-mass energy. The small 
variations of the beam energy over the duration of the run 
are corrected when calculating toes- For the momentum 
Pi {i = 0,B) and the energy Eq, the subscripts and B 
refer to the e~'"e~ system and the reconstructed B meson, 
respectively. The energies -B^(.) and E^ are calculated 
from the measured £)(*)" and h° momenta. 

For the various decay channels of the B signal events, 
the Toes distribution peaks at the B mass with a resolu- 
tion of 2.6 — 3 MeV/c^, dominated by the beam energy 
spread, whereas AE peaks near zero with a resolution of 
15 — 50 MeV depending on the number of photons in the 
final state. 



The _D*" mesons are reconstructed in the I?°7r° and 
D^j decay modes. The tt" and D" candidates arc re- 
quested to satisfy the selections described in Sec. IIIIB II 
and nil B7l respectivelv. The photons from D*'^ — > D'^7 
must fulfill the additional condition £^7 > 130 MeV and 
must pass the tt" veto as described in Sec. IIIIB 21 

The resolution of the mass difference Am = mjyto — 
rriDO is about 1.3 MeV/c^ for D^n" and 7 MeV/c^ for 



C. Selection of B-meson candidates 

The B candidates are reconstructed by combining a 
D^*)" with an /i°, with the D'^*^° and h° masses con- 
strained to their nominal values (except when h^ is an 
to). One needs to discriminate between true B signal 
candidates and fake B candidates. The fake B candi- 
dates originate from combinatorial backgrounds, from 
other specific B modes, or from the cross feed events 
between reconstructed color-suppressed signals. 



1. B -mesons kinematic variables 

Two kinematic variables are commonly used in BaBAR 
to select B candidates: the energy-substituted mass toes 
and the energy difference AE. These two variables use 
constraints from the precise knowledge of the beam ener- 
gies and from energy conservation in the two-body decay 
T(45) — ?> BB. The quantity tties is the invariant mass of 



2. Rejection ofe'^t 



qq background 



The continuum background e^e^ — > qq, where q is a 
light quark m, d, s, or c, creates high momentum mesons 
£)(*)o^ 7r°, ri^\ oj that can fake the signal mesons origi- 
nating from the two body decays B° — > D^*^°h^. That 
background is dominated by cc processes and to a lesser 
extent by ss processes. Since the B mesons are produced 
almost at rest in the T(4S') frame, the r(4S') -^ BB 
event shape is isotropically distributed. By comparison, 
the qq events have a back-to-back jet-like shape. The 
qq background is therefore discriminated by employing 
event shape variables. The following set of variables was 
found to be optimal among various tested configurations: 

• The thrust angle 9t defined as the angle between 
the thrust axis of the B candidate and the thrust 
axis of the rest of event, the thrust axis being the 
axis on which the sum of projected momentum is 
maximal. The distribution of | cos 9t \ is flat for 
signal and peaks at 1 for continuum background. 



Event shape monomials Lq and L2 defined as 



Lo = Y.\-f*\; L2 = Y.\f* 



(3) 



with p* being the CM momentum of the particle 
i that does not come from the B candidate, and 9* 
is the angle between p * and the thrust axis of the 
B candidate. 
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• The polar angle 9g between the B momentum in 
the T(45) frame and the beam axis. With the 
T{AS) being vector and the B mesons being pseu- 
doscalar, the angular distribution is proportional to 
sin^ 9g for signal and roughly flat for baekground. 

These four variables are combined into a Fisher dis- 
criminant built with the TMVA [3J| toolkit package. An 
alternate approach employing a multi-layer perceptron 
artificial neural network with two hidden layers within 
the same framework was tested and showed marginal rel- 
ative gain, therefore the Fisher discriminant is used. 

The Fisher discriminant J-"shapc is trained with signal 
MC events and off-peak data events. In order to maxi- 
mize the number of off-peak events all the B^ — >■ D^*^'^hP 
modes are combined. We retain signal MC events with 
TTiES in the signal region 5.27—5.29 MeV/c^ and off-peak 
data events with tties in the range 5.25 — 5.27 MeV/c^, 
accounting for half of the 40 MeV CM energy-shift be- 
low the T{AS) resonance. The training and testing 
of the multivariate classifier are performed with non- 
overlapping data samples of equal size obtained from a 
cocktail of 20, 000 MC simulation signal events and from 
20, 000 off-peak events. The obtained Fisher formula is 



J-^hapc = 2.36 - 1.18 X |cos6't| + 
0.20 X io - 1.01 X La - 0.80 x | cos^lj] 



(4) 



The qq background is reduced by applying a selection cut 
on J- shape- The selection is optimized for each of the 72 
signal decay channels by maximizing the statistical sig- 
nificance with signal MC against generic MC e+e" — >■ qq, 
q =/: b. This requirement for the various decay modes re- 
tains between about 30% and 97% of B signal events, 
while rejecting between about 98% and 35% of the back- 
ground from light qq. 



3. Rejection of other specific backgrounds 



We require B° — >■ D^u candidates to satisfy J^hci > —0.1, 
to obtain an efficiency (rejection) on signal (background) 
of about 85% (62%). 

We also exploit the angular distribution properties in 
the decay D*" — > D^tt" to reject combinatorial back- 
ground. We define the helicity angle 9]j* as the angle 
between the line-of-fiight of the -D° and that of the i?°, 
both evaluated in the D*° rest frame. The angular distri- 
bution is proportional to eos^ 9d* for signal and roughly 
flat for combinatorial background. Although in princi- 
ple such a behavior could be employed for B'^ -> £'*°7r°, 
D*^r], and D*'^ij', a selection on |cos6'£). | significantly 
improves the statistical significance for the B^ — > £'*°7r° 
mode only. Therefore D*° candidates coming from the 
decay B^ — )• D*°7r° arc required to satisfy | cos6'£i. | > 0.4 
with an efficiency (rejection) on signal (background) of 
about 91% (33%). 

A major BB background contribution in the analy- 
sis of the B^ — > £)''*-*°7r'' decay channel comes from the 
color-allowed decay B~ — >■ D^*^^p~ . If the charged pion 
(mostly slow) from the decay p^ — > 7r~7r° is omitted in 
the reconstruction of the B° candidate, B^ -^ D^*^^ p'' 
events can mimic the L)(*)°7r'^ signal. Moreover, the de- 
cay modes B{B^ — > Z3(*)°p^) are 30 — 50 times larger 
than those of the B^ — > D(*'°7r° modes, and are poorly 
known: AB/B = 13.4%- 17.3% [23_ A veto is applied to 
reduce this background. For each B'^ — >■ r)(*)°7r° candi- 
date, we combine any remaining negatively charged track 
in the event to reconstruct a B~ candidate in the decay 
mode D'^*^^ p" . If the reconstructed B~ candidate sat- 
isfies mEs(S-) > 5.27 GeV/c2, |A£;(E-)| < 100 MeV, 
and \mp- - m^P^I < 250 MeV/c^, then the initial B° 
candidate is rejected. For the analysis of the decay mode 
^° -^ D°TT° {B- -J> D*°Tr°), the veto retains about 
90% (82%) of signal and rejects about 67% (56%) of 



B- ^ D^p- and 44% (66%) of B' 
ground. 



D* 



back- 



The Lu mesons in i?° — )• D'^uj decays arc longitudinally 
polarized. We define the angle 9i^ [3, [SJI as the the angle 
between the normal to the plane of the three daughter 
pions in the uj frame and the line-of- flight of the B'^ meson 
in the uj rest frame. This definition is the equivalent 
of the two-body helicity angle for the three-body decay. 
To describe the three-body decay distribution of w — > 
TT'^TT~TT^, we define the Dalitz angle 9d [7| as the angle 
between the 7r° momentum in the a; frame and the 7r+ 
momentum in the frame of the pair of charged pions. 

The signal distribution is proportional to cos^ 9^^ and 
sin 9o, while the combinatorial background distribution 
is roughly flat as a function of cos 9^ and cos 9]j . These 
two angles are combined into a Fisher discriminant J-hd 
built from signal MC events and generic qq and BB MC 
events: 



j:^^!^ -1.41-1.01 X |cos6i£,| +3.03 x |cos( 



(5) 



4. Choice of the "best" B candidate in the event 

The average number of B'^ — > D'-*^'^h^ candidates per 
event after all selections ranges between 1 and 1.6 de- 
pending on the complexity of the sub-decays. We per- 
form all the 72 B^ — >■ D'^*^^h^ decay mode analyzes in 
parallel. Such that each possible decay channel is se- 
lected with a dedicated analysis, for a given decay mode 
one B candidate only is kept per event. The chosen B is 
that with the smallest value of 



^^ _ / mpo - fim^o \ ' ^ / m,,o - ^im,„ \ ^ ^ ^^^ 



a^ 



rrioO 



V ^'■"1^0 
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for D'^h^ modes, and 



Xb 



2 _ frriDO - fl,n^o \ ^ /™/i« - A^™,.o ^ ^ 



/Am — /lA ^ 

\ CTAm 



'™hO 



(7) 



for the D*^hP modes. The quantities cfm o ^'^d Um o 
(/i„i^o ^nd /im^o ) ^'''6 the resolution (mean) of the mass 
distributions. The quantities ^Am and CTAm are respec- 
tively the mean and resolution of the Am distributions. 
These quantities are obtained from fits of the mass distri- 
bution of simulated candidates selected from signal MC 
simulations. 

The probability of choosing the true B^ candidate in 
the event according to the above criteria ranges from 71 
to 100%. The cases with lower probabilities correspond 
to the D(*)°/i° modes with high neutral multiplicity. 



5. Selection efficiencies 

The branching fractions of the B^ — > D^*^^hP decays 
is computed as 



6(5° ^ i?(*)"/i°) 



Ns 



Nbb-£- Bs. 



(8) 



where Bscc is the product of the branching fractions as- 
sociated with the secondary decays of the Z)'^*^'^ and /i° 
mesons for the each of the 72 decay channels considered 
in this paper [2g|. N^g is the number of BB pairs in 
data and Ns is the number of signal events remaining 
after all the selections. The quantity £ is the total signal 
efhciency including reconstruction (detector and trigger 
acceptance) and analysis selections. It is computed from 
each of the 72 exclusive high statistics MC simulation 
samples. 

The selection efficiency from MC simulation is slightly 
different from the efficiency in data. The MC efficiency 
and its systematic uncertainty therefore has to be ad- 
justed according to control samples. For the reconstruc- 
tion of 7r°/7, the efficiency corrections are obtained from 
detailed studies performed with a high statistics and 
high purity control sample of 7r° mesons produced in 
T — >■ p{t:t:'^)vt decays normalized to r — >■ i:Vt, to unfold 
tracking effects. Such corrections are validated against 
studies performed on the relative ratio of the number of 
detected D^ mesons in the decays D° -^ K~tt~^7t'^ and 
Z?" -^ Jv ^7r+, and produced in the decay of £)*+ mesons 
from e^e^ -^ cc events. The relative data/simulation 
efficiency measurements for charged tracks are simi- 
larly based on studies of track mis-reconstruction us- 
ing e'^e~ — > T^r~ events. On one side the events are 
tagged from a lepton in the decay t~ — >■ 1~uiVt and on 
the other side one reconstructs two or three tracks from 
the decay r"*" — >■ ■K'^T:~h^VT. The simulated efficiency 



of charged particle identification is compared to the ef- 
ficiency computed in data with control samples of kaons 



from D*+ -^ £)°(iir^ 7r+)7r+ produced in e+e^ — > cc 
events. The efficiency for Xg candidates is modified using 
a data sample of K^, mainly arising from the continuum 
processes e^ e" into qq. 

The efficiency corrections for the selection criteria ap- 
plied to D^*'^ candidates and on the Fisher discriminant 
(-Tiici) for the continuum qq {q ^ b) rejection are obtained 
from studies of a B~ — >■ £>(*^°7r~ control sample. This 
abundant control sample is chosen for its kinematic simi- 
larity with B^ — > D^*^^h^. The corrections are computed 
from the ratios £'rei.(data)/frei.(MC), where the relative 
efficiencies frei. are computed with the signal yields as ob- 
tained from fits to rrtEs distributions of B~ — >■ D^*^^t:~ 
candidates in data and MC simulation, before and af- 
ter applying the various selections. The obtained re- 
sults are checked with the color-allowed control sample 
B^ — > D^*^^p^, which has slightly different kinematics 
due to the relatively higher mass of the p~ , and therefore 
validates those corrections for the modes such as D'^*^'^rj' . 

The reconstruction efficiency of B^ — > D*^ui depends 
on the angular distribution, which is not yet known 
(/l ~ 0.5 — 1). To evaluate this efficiency we combine 
a set of properly weighted fully longitudinally and fully 
transversely polarized MC samples, according to the frac- 
tion of longitudinal polarization that we measure in this 
paper (see Sec. W^ . 



D. Fit procedure and data distributions 

We present the fits used to extract the branching frac- 
tions B. For each of the 72 possible B° -^ D'-*^^h° sub- 
decay modes, using an iterative procedure (discussed in 
Sec. IIIID5|) . we fit the AE distribution in the range 
\AE\ < 280 MeV for rngs > 5.27 GeV/c^ to get the 
signal (Ng) and background yields. The fit of the AE 
distribution allows us to model and adjust the complex 
non-combinatoric B background structure without rely- 
ing completely on simulation. 

The data samples corresponding to each B^ decay 
mode are disjoint and the fits are performed indepen- 
dently for each mode. According to their physical ori- 
gin, four categories of events with differently shaped AE 
distributions are considered: signal events, cross-feed 
events, peaking background events, and combinatorial 
background events. The event (signal and background) 
yields are obtained from unbinned extended maximum 
likelihood (ML) fits. We write the extended likelihood £ 
as 



N 



L = 



-,N 



N\ 



n/(Ai?ji^,"), 



(9) 



where 9 indicates the set of parameters which are fitted 
from the data. TV is the total number of signal and back- 
ground events for each sub-decay mode, and n = ^^ Ni is 
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FIG. 2: Fits of A£; distributions in data for modes B° -^ D%° (a), B° -^ D°uj (b), B° -^ D°ri{'y-f) (c), B° -^ D°Ti{TTmr°) (d), 
B" — >■ D°ri'{TVTvrj) (e), and B" — > D^rj^p^j) (f). The data points with error bars are measurements in data, the curves are the 
various PDF components: the soUd (blue) fitted total PDF, the dotted (red) signal PDF, the dotted-dashed (black) cross-feed 
PDF, the double dotted-dashed (brown) B~ — >■ D^*'^'p~ PDF, and the long dashed (blue) combinatorial background PDF. 



the expectation value for the total number of events. The 
sum runs over the different expected number Ni of signal 
and background events in the various i categories. The 
total probability density function (PDF) f{AEj\6,n) is 
written as the sum over the different signal and back- 
ground categories, 



fiAE,\e,n) 



j:^N,MAE,\e) 



(10) 



where fi{AE\9) is the PDF of category i (signal or back- 
ground component). Some of the PDF component pa- 
rameters are fixed from the MC simulation (see details 
in the following sections). 

The individual corresponding branching ratios are 
computed and then combined as explained in Sec. |Vl 



1. Signal contribution 

All of the 72 possible reconstructed i?° decay channels 
contain at least one photon. Due to the possible energy 
losses of early showering photons in the detector mate- 
rial before the EMC, the AE shape for signal is modeled 
by the modified Novosibirsk PDF |3l|. A Gaussian PDF 
is added to the modes with a large AE resolution to 
describe mis-reconstructed events. The signal shape pa- 
rameters are estimated from a ML fit to the distributions 



of simulated signal events in the high statistics exclusive 
decay modes. 



2. Cross-feed contribution 

We call "cross feed" the events from all of the recon- 
structed D^*^^h^ modes, except the one under consider- 
ation, that pass the complete selection. The cross-feed 
events are a non-negligible part of the AE peak in some 
of the modes, and the signal event yield must be corrected 
for these cross-feed events. As the various decay channels 
are studied in parallel, we use an iterative procedure to 
account for those contributions in the synchronous mea- 
surements (see Sec. IIIIETS]) . 

The dominant cross- feed contribution to 5° -^ D'^h^ 
comes from the companion decay channel B" -^ D*^hP ^ 
when the 7r°/7 from the D*^ decay is not reconstructed. 
Such cross-feed events are shifted in AE by approxi- 
mately the mass of the tt" (—135 MeV), with a long 
tail from Z?*°(— >■ D°'y)h'^ leaking into the signal region. 
Similarly, the decay channel B" — > D*'^h^ receives a 
cross-feed contribution from the associated decay mode 
B^ -^ D^hP and there is a cross-contamination between 
the D*° -^ D°Tr° and D*° -^ D°j decay channels. 

The remaining cross- feed contributions, i.e. from other 
B^ -^ D'^*^^hP color-suppressed decay modes, were stud- 
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FIG. 3: Fits of AE distributions in data for modes S" — > D* 



B" -^ D'^^u (b), B° -^ I>*°77(77) (c), B° 



D'^riimv-K") 



(d), B° -s- D*°rj'(nnri) (e), and B" -^ D'^rj'ip^j) (f), where the D*" mesons decay into the signal mode D°n°. A detailed 
legend is provided in the caption of Fig. |21 



led with the generic MC simulation. They were found to 
be highly negligible and represent at most 1 % of the sig- 
nal, in the region \AE\ < 100 MeV. Therefore they are 
accounted for by the generic MC simulation, whose gen- 
erated branching fractions were taken from the PDG [23 . 



3. Peaking BE background contributions 



The major background in the reconstruction of S" — >■ 
J5(*)07^o comes from the decays B~ — >■ D'^*^"p~ 
(see Sec. IIIIC3p . Their contribution is modeled by 
histogram-based PDFs built from the high statistics ex- 
clusive signal MC simulation samples. The individ- 
ual distributions of the two backgrounds B^ — > D^ p^ 
and B- -^ D*°p- that pass the B° -)> D<^*^°tt° selec- 
tions, including the specific veto requirement as described 
in Sec. IIII C 3[ cannot be distinguished. As a conse- 
quence, given the large uncertainty on their branching 
fractions, the overall normalization of B~ —J- D^*^^ p~ 
PDF is left floating but the relative ratio N{B^ ^■ 
D*°p-)/N{B- -^ D^p~) of the PDF normalization 
is fixed. The value of this ratio is extracted directly 
from the data by reconstructing exclusively each of the 
B~ — 7> D'^*^^ p~ modes rejected by the veto require- 
ments. Those fully reconstructed B~ mesons differ from 
the B- ^ D^*^°p-, that pass all the B° -^ D^^'Ott" 
selections, by the additional selected slow charged tt 



originated from the p~ meson. The relative correction 
on that ratio for events surviving the veto selection is 
then computed using the MC simulation for generated 
B^ — > D^*^^p^ decays. A systematic uncertainty on that 
assumption is assigned (see Sec. |TV|. 

In the cases of B° — >• D(*)°a;/77(— > 7r+7r^7r°) modes, 
additional contributions come from the B decay modes 
D^*'mnT^'^' , where n = 1, 2, or 3, and through interme- 
diate resonances such as to and pg (1690)(— > ujn^). These 
peaking backgrounds are modeled by a first-order poly- 
nomial PDF plus a Gaussian PDF determined from the 
generic BB MC simulation. The relative normalization 
of that Gaussian PDF component is left floating in the 
fit, since some of the branching fractions of the B decay 
modes D^*") nnn^^^ are not precisely known |2a |. 



4-. Combinatorial background contribution 

The shape parameters of the combinatorial back- 
ground PDFs are obtained from ML fits to the generic 
BB and continuum MC, where all signal, cross feed 
and peaking BB background events have been removed. 
The combinatorial background from BB and qq {q 7^ h) 
are summed and modeled by a second-order polynomial 
PDF. 
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FIG. 4: Fits of AE distributions in data for modes B° -^ D'°tv° (a), B° -^ D'^oj (b), B° -^ D''°-n{-yy) (c), B° -^ D*°rj(Tv-Kn°) 
(d), and B" — >■ D*° rj' {-rnvrj) (e), where the D*" mesons decay into the signal mode D^^j. The unfitted AE distribution of 
B" — >■ D*^ [D'^ ^)ri' {p° ')) candidates is also displayed (f). A detailed legend is provided in the caption of Fig. [21 



5. Iterative fitting procedure 

We fit the /S.E distribution using the PDFs for the sig- 
nal, for the cross feed, for the peaking background, and 
for the combinatorial background as detailed in the pre- 
vious sections. The normalization of the signal, the peak- 
ing BB backgrounds, and the combinatorial background 
components arc allowed to float in the fit. The mean of 
the signal PDF is left floating for the sum of £)'*'" sub- 
decays. For each D° sub-mode, the signal mean PDF 
is fixed to the value obtained from the fit to the sum of 
D^ sub- modes. Those free parameters are extracted by 
maximizing the unbinned extended likelihood to the A_E 
distribution defined in Eqs. ^ and ^. Other PDF 
parameters are fixed from fit results obtained with MC 
simulations, when studying separately each of the signal 
and background categories. 

In the global event yield extraction of all the various 
i3° — >■ Z?(*)°/i° color-suppressed signals studied in this 
paper, a given mode can be signal and cross feed to other 
modes at the same time. In order to use the B computed 
in this analysis, the yield extraction is performed through 
an iterative fit on D*'^h^ and D^h^ . The normalization 
of cross-feed contribution from D^*^°h^ is then fixed to 
the B measured in the previous fit iteration. For the 
cross-feed contributions, the PDG branching fraction [2Q| 
values are used as starting points. This iterative method 
converges quickly to stable B values, with a variation of 



less than 10% of the statistical uncertainty, in less than 
5 iterations. 

We check the absence of biases in our fit procedure 
by studying pseudo-experiments with a large number of 
different data-sized samples for the various signals. The 
extraction procedure is applied to these samples where 
background events are generated and added from the fit- 
ted PDFs. The signal samples are assembled from non- 
overlapping samples corresponding to the exclusive high 
statistics MC signals, with yields corresponding to the 
MC-generated value of the branching fraction. No signif- 
icant biases are found. 



6. Data distributions and event yields from summed 
sub-decay modes 

The fitting procedure is applied to data at the very 
last stage of the blind analysis. Though the event yields 
and B measurements are performed separately for each 
of the 72 considered sub-decay modes, we illustrate here, 
in a compact manner, the magnitude of the signal and 
background component yields, and the statistical signif- 
icances, of the various decay channels -B" -^ D^*'>'^h'^, 
summing together all the D° sub-modes. The fitted AE 
distributions, for the sum of D° sub-modes, are given 
in Figs, mil andH for, respectively, the B° -^ D^h", 
D*°{^ D"n")h°, and D*°{-^ D^j)h^ modes. 
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TABLE I: Numbers of signal events (Ns), combinatorial background (Ai'combi), cross feed (Net), and B~ — >■ D^^'^ p~ {N-op) 
events computed from the AE fits to data and counted in a signal box \AE\ < 2.5 a, together with the statistical significances 
in numbers of standard deviations 5stat (see text). The quoted uncertainties are statistical only. 



B° ^ 


Ns 




-^^coinb 




A^cf 




Nop 


Statistical 


(decay channel) 
























significance 5stat 


D^o 


3429 


± 


123 


2625 


± 


75 


97 


± 


3 


700 


± 14 


41 


D^vi-yy) 


1022 


± 


55 


532 


± 


14 


13 


± 


1 




- 


36 


D°r,(-KTVK°) 


411 


± 


29 


191 


± 


6 


2 


± 







- 


23 


DOlo 


1374 


± 


120 


886 


± 


25 


18 


± 


2 




- 


38 


D^n' {ir-Krii-n)) 


122 


± 


13 


41 


± 


3 




- 






- 


14 


D%'{p°j) 


234 


± 


40 


1253 


± 


17 


1 


± 







- 


7.4 


D*0(D«7rO)7rO 


883 


± 


40 


268 


± 


21 


39 


± 


2 


175 


± 5 


34 


D*0{DO'y)nO 


622 


± 


47 


469 


± 


33 


295 


± 


23 


602 


± 20 


17 


D*0(Z)0^0)r,(77) 


338 


± 


25 


201 


± 


9 


17 


± 


1 




- 


19 


D*0(Z)07)r,(77) 


187 


± 


24 


254 


± 


12 


85 


± 


11 




- 


8.7 


D*0[D°TrO)n{nnnO) 


123 


± 


15 


90 


± 


4 


5 


± 


1 




- 


11 


D*0{DO'y)r,{nW) 


88 


± 


14 


65 


± 


4 


16 


± 


3 




- 


7.6 


D'0{D'>nO)LU 


806 


± 


48 


1365 


± 


18 


33 


± 


2 




- 


20 


D''0{D°'y)LU 


414 


± 


44 


1290 


± 


19 


132 


± 


14 




- 


10 


D*0(D0 7rO)r;'{7r7r7?) 


45 


± 


8 


18 


± 


2 


2 


± 







- 


8.5 


D*"{D"^W{mTri) 


12 


± 


5 


8 


± 


1 


5 


± 


2 




- 


3.2 


D*0{DOnO)r,'{pO-y) 


115 


± 


25 


487 


± 


11 


3 


± 


1 




- 


5.4 



The signal and background yields obtained from the 
fit to the summed sub-mode data for the 5° — J> D^*^°h° 
are presented in Table U with the corresponding statis- 
tical significances. The signal and background yields are 
computed in the signal region \AE\ < 2.5 a (where a is 
the signal resolution). In the same range we calculate 
the statistical significance of the various signals from the 
cumulative Poisson probability p to have a background 
statistical fluctuation reaching the observed data yield, 



-t-oo 

E 

fc=Wca„d 



fc! 



(11) 



where A'cand is the total number of selected candidates 
in the signal region and v the mean value of the total 
expected background, as extracted from the fit. This 
probability is then converted into a number of equivalent 
one-sided standard deviations Sstut, 



Sstut = V2 erf cInverse(p/2). 



(12) 



The function erf cinverse is the inverse of the comple- 
mentary error function (see statistics review in j2a|). 

The majority of the decay channels present clear 
and significant signals. In particular, the modes 
ZJ^Ty' (777777(77)) ^-nd D^ri'{p^j) arc observed for the first 
time. 

Before performing the final unblinded fits on data, 
among the various 72 initial possible decay channels, sev- 
eral sub-decay modes have been discarded. The decision 
to remove those sub-modes have been taken according 
to analyses performed on MC simulation, as no signifi- 
cant signals are expected. The ehminated decay channels 



are: B° -^ D^*^°r]' and D*^{D^j)T]{Tnnr°), where D° -^ 
K°n+n-, D*°{D°-f)r]'{Tnnj), where D° -^ K--k+it-t^+, 
as well as the whole decay channel D*^ [D^ ^)rj' [p^ ^) . 
These are sub-modes with poor signal efficiency, caused 
by large track multiplicity or modest D^ secondary 
branching fractions, such that the expected signal yields 
are very low. In addition, they have large background 
contributions. We concluded that adding such decay 
channels in the global combinations would degrade the 
B measurements. These choices based on a Monte Carlo 
simulation only studies have been confirmed in data (see 
for example Fig. U (bottom right)). 



IV. SYSTEMATIC UNCERTAINTIES ON 
BRANCHING FRACTIONS 

There are several possible sources of systematic uncer- 
tainties in this analysis, which are summarized in Ta- 
bic HI 

The categories "77^/7 detection" and "Tracking" ac- 
count respectively for the systematic uncertainties on the 
reconstruction of 77*^/7 and for charged particle tracks, 
and are taken from the efficiency corrections computed 
in the studies of r decays from e^e~ — >■ t'^t~ events (see 
Sec. lIIICSI) . 

Similarly, the systematic uncertainties on kaon iden- 
tification and on the reconstruction of K^ mesons are 
estimated from the MC efficiency corrections computed 
in the study of pure samples of kaons and K^ mesons 
compared to data (see Sec. IIII(j5l) . 

The uncertainty on the secondary branching fraction 
results from our limited knowledge of the D^*'^ and h^ 
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TABLE II: Combined contributions to the branching fraction S(-B" — )> Z?'*^"/!") relative systematic uncertainties (%). 



Sources 








AB/B{%) for the B' 


' decay 












D^tt" 


D°vill) 


DO-niTTTTTrO) 


DOl, 


D°ri'{TTnri) 


D' 


'V{pS) 


D'Ott" 


D-^'r, 


D'^'lu 


D*%' 


■rfi 1^ detection 


3.5 


3.5 


3.6 


3.6 


3.7 






2.3 


6.2 


5.5 


5.7 


5.8 


Tracking 


0.9 


0.9 


1.6 


1.7 


1.6 






1.6 


0.9 


1.1 


1.6 


1.6 


Kaon ID 


1.0 


1.1 


1.1 


1.1 


1.2 






1.1 


1.1 


1.1 


1.1 


1.2 


K% reconstruction 


0.7 


0.7 


0.6 


0.8 


0.6 






0.6 


0.4 


0.3 


0.5 


- 


Secondary B 


1.6 


1.6 


2.0 


1.8 


2.3 






2.4 


5.1 


5.7 


5.5 


5.1 


BB counting 


1.1 


1.1 


1.1 


1.1 


1.1 






1.1 


1.1 


1.1 


1.1 


1.1 


MC statistics 


0.1 


0.2 


0.3 


0.2 


0.4 






0.4 


0.2 


0.2 


0.3 


0.3 


Particles selection 


0.3 


0.4 


0.2 


1.0 


0.3 






1.0 


0.2 


0.1 


0.1 


1.2 


AEfit 


2.1 


2.2 


1.3 


2.1 


1.1 






2.1 


1.0 


0.6 


1.4 


0.5 


Combinatorial background 


0.1 


0.1 


0.1 


0.1 


0.1 






0.1 


0.1 


0.1 


0.1 


0.1 


D(*)°p- background 


1.8 


- 


- 


- 


- 






- 


5.6 


- 


- 


- 


D*^u} polarization 


- 


- 


- 


- 


- 






- 


- 


- 


1.4 


- 


Total 


5.1 


4.9 


4.9 


5.3 


5.1 






4.7 


9.6 


8.2 


8.5 


8.2 



TABLE III: Branching fractions of the decay channels B" — )> D^*'"hP measured in the different secondary decay modes. The 
first uncertainty is statistical and the second is systematic. The cells with "-" correspond to decay channels that have been 
discarded after the analysis on simulation, and confirmed with data, as no significant signal is expected or seen for them. 



B(BO^) (xlO-4) 


D" ^ Kit 


D" -> K3-7T 


DO ^ Ktttt" 


D" -^ KOtt+tt- 


DOttO 


2.49 ± 0.13 ± 0.16 


2.69 ± 0.15 ± 0.17 


2.97 ± 0.15 ± 0.25 


2.90 ± 0.28 ± 0.23 


D0r,(77) 


2.46 ± 0.18 ± 0.14 


2.56 ± 0.19 ± 0.16 


2.37 ± 0.20 ± 0.20 


2.62 ± 0.37 ± 0.21 


DO»7(7r7r7rO) 


2.59 ± 0.27 ± 0.12 


2.65 ± 0.30 ± 0.14 


2.48 ± 0.29 ± 0.20 


2.28 ± 0.54 ± 0.18 


D^UJ 


2.59 ± 0.18 ± 0.20 


2.34 ± 0.19 ± 0.15 


2.42 ± 0.20 ± 0.21 


3.17 ± 0.39 ± 0.24 


DOn'iniTTji'yy)) 


1.40 ± 0.25 ± 0.07 


1.37 ± 0.26 ± 0.08 


1.34 ± 0.27 ± 0.11 


1.30 ± 0.50 ± 0.12 


D^v'ip"-/) 


1.58 ± 0.42 ± 0.09 


1.79 ± 0.57 ± 0.10 


1.91 ± 0.54 ± 0.15 


1.55 ± 0.89 ± 0.16 


D*'\D°nO)n'> 


2.95 ± 0.25 ± 0.30 


2.95 ± 0.29 ± 0.33 


3.52 ± 0.29 ± 0.43 


2.32 ± 0.56 ± 0.24 


D*0(D07)t° 


3.49 ± 0.40 ± 0.83 


2.25 ± 0.50 ± 0.63 


3.02 ± 0.50 ± 0.90 


3.53 ± 1.14 ± 0.99 


D'-0{D°7T'>)r,{-yy) 


2.52 ± 0.32 ± 0.26 


2.57 ± 0.33 ± 0.29 


2.41 ± 0.32 ± 0.32 


4.09 ± 0.74 ± 0.49 


D*0(DS)r,(77) 


2.62 ± 0.45 ± 0.33 


2.81 ± 0.49 ± 0.35 


2.87 ± 0.55 ± 0.39 


2.75 ± 0.78 ± 0.36 


D*0{DOTTO)r,{mTn°) 


2.27 ± 0.50 ± 0.20 


2.60 ± 0.55 ± 0.24 


1.93 ± 0.46 ± 0.22 


1.21 ± 0.87 ± 0.13 


D'^{D°-i)r,(-KTVK^) 


2.93 ± 0.71 ± 0.32 


2.55 ± 0.80 ± 0.29 


1.94 ± 0.81 ± 0.24 


- 


D'0{DOn'>)iU 


5.07 ± 0.45 ± 0.47 


4.00 ± 0.49 ± 0.36 


4.38 ± 0.51 ± 0.51 


5.02 ± 0.98 ± 0.53 


D*^{D°-t)ui 


3.66 ± 0.64 ± 0.41 


4.46 ± 0.80 ± 0.56 


4.59 ± 0.87 ± 0.57 


4.28 ± 1.71 ± 0.57 


D'0{DOnOW{nnv{^^)) 


1.09 ± 0.38 ± 0.09 


1.67 ± 0.44 ± 0.15 


1.34 ± 0.49 ± 0.15 


- 


D*0(D07)r,'(7r7rr,(77)) 


0.75 ± 0.49 ± 0.24 


- 


1.19 ± 0.69 ± 0.39 


- 


D*0(D07rO)r?'(p07) 


2.10 ± 0.82 ± 0.23 


1.21 ± 0.90 ± 0.14 


1.45 ± 0.95 ± 0.18 


- 


D-0(D07)r,'(p07) 


- 


- 


- 


- 



sub-mode branching ratios [26| (including secondary de- 
cays into detected stable particles) . Correlations between 
the different decay channels were accounted for. 

The uncertainty related to the number of BB pairs and 
the limited available MC-sample statistics when comput- 
ing the efficiency of various selection criteria are also in- 
cluded. 

Systematic uncertainties due to the intermediate par- 
ticles mass selections are computed as the relative differ- 
ence of signal yield when the values of the mass means 
and mass resolutions are taken from a fit to the data. 
Systematic effects from the qq {q ^ b) rejection and 
the Z?(*'° selections are obtained from the study per- 



formed on the control sample B~ — !• D^*^°7r~ and are 
estimated as the limited confidence on the efficiency cor- 
rection ratio: £roi.(data)/£i.ci.(MC), including the corre- 
lations between the samples before and after selections 
(see Sec. IIIIC Sp . The effects of the cuts on p" and D'^oj 
helicities are obtained by varying the selection cut values 
by ±10% around the maximum of statistical significance. 
All uncertainties on intermediate particle selections are 
combined into the category "Particles selection" . 

The uncertainty quoted for "AE Fit" gathers the 
changes due the limited knowledge on the shapes of signal 
and background PDFs, and on the cross-feed branching 
fraction. For example, the modes D'^*^"7r°/r/(77) have 
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high or low momcntuni 7 in the final state. The difference 
between data and MC simulation in energy scale and res- 
olution for neutrals is estimated from a study of the high 
statistics control sample B~ — )- D^{K~7T~^)p~ {tt^tt~), 
which yields the difference between data and MC sim- 
ulation, ~ 5.7 MeV, for the mean and ~ 3.3 MeV, for 
the resolution. This study was cross-checked against an- 
other one performed with the high statistics and very 
pure control sample B~ -^ D^ {K~ tt^ tt^)-!!' . With that 
control sample we find a difference between data and 
MC simulation of respectively ~ 1.7 MeV, for the mean 
and ~ 0.04 MeV, for the resolution. For the modes 
B^ — >■ I?(*)°7r''/77(77), the uncertainty due to the signal 
shape is obtained by varying conservatively the signal 
PDF mean by ±5.7 MeV and the width by ±3.3 MeV. 
For the other B^ signal modes, each PDF parameter is 
varied within ±1 cr of its MC simulation precision, and 
the relative difference on the fitted event yield is taken as 
a systematic incertitude. The various parameters are var- 
ied one at a time, independently. The relative differences 
while varying the IS.E PDF parameters are then summed 
up in quadrature. This sum is taken as the systematic 
doubt on the /S.E shapes. 

The uncertainty on the combinatorial background 
shape (including both BB and qq {q 7^ b) events) is eval- 
uated from the comparison of generic MC simulation and 
data in the rriES sidebands: 5.24 < tties < 5.26 GeV/c^. 
The observed difference between data and simulation has 
then been used as a systematics. When a Gaussian is 
added to the combinatorial background shape, to model 
additional peaking BB background contributions (see 
Sec. IIIID3|) . the related effect is computed by varying 
its means and resolution by ±1 a. 

We account for possible differences in the PDF shape 
of the B~ — >■ D^*^'^p^ background that is modeled by a 
non-parametric PDF. As above, it is obtained by shifting 
and smearing the PDF mean and resolution by ±5.7 MeV 
and ±3.3 MeV respectively. The non-parametric PDF is 
convoluted with a Gaussian with the previously defined 
mean and width values. The quadratic sum of the var- 
ious changes on the signal event yield is taken as the 
systematic uncertainty. 

The relative ratio of the B~ — >• D*°p~ and B~ — ^ 
D^ p~ backgrounds for the studies of the modes B^ —>■ 
D(*)"7r" has been fixed to the selected B- -^ D'^*^°p- 
events of the data control sample, for rejected B^ events 
with the veto described in Sec. IIII C 31 The effect of such 
a veto on that ratio is then computed from MC simula- 
tion. We assign as a conservative systematic uncertainty 
half of the difference between the nominal result and 
the result from the MC simulation assuming the PDG 
branching ratios of B' -> D'^*^°p- ^. 

The acceptance of i?° -> I?*°Cti is computed from the 
sum of purely longitudinally (/l = 0) and transversely 
(/l = 1) polarized MC simulation signals, weighted by 
our measurement of /^ (see Sec. lVI|) . The systematic lim- 
ited knowledge of the efficiency due to the unknown frac- 
tion of D*'^Lu longitudinal polarization is then estimated 



by varying /^ by ±1 cr in the estimation of the signal 
acceptance. This contribution is slightly larger than 1%, 
while it is expected to be about 10.5% if the fraction /l 
were unknown. This is one of the motivations for mea- 
suring the polarization of the decay channel B" — >■ £)*°w 

(sec Sec. El)- 

The most significant sources of systematic uncertain- 
ties come from the 7r*'/7 reconstruction, the AE fits, 
and the uncertainties on the world average branching 
fractions of the secondary decay channels. In the case 
of the modes B^ — >■ D(*)°7r°, the contributions from 
B^ -^ D^*^^p~ backgrounds arc also not negligible. 



TABLE IV: Branching fractions of decay channels B — >■ 
£j{»)o^o^ where the branching fraction measured in each D° 
modes are combined. For the modes with h = 77, r;', we give 
the combination (comb.) of the branching fraction computed 
with each sub-modes of 77' ' . Tire first uncertainty is statistical 
and the second is systematics. The quality of the combination 
is given through tire value of x^/ndof , with tire corresponding 
probability (p-value) given in parenthesis in percents. 



B° mode 


Z3(xlO-^) 


X^/ndof 
(p-value %) 


i3"7r" 


2.69 ± 0.09 ± 0.13 


2.81/3 (42.2) 


D°r} (comb.) 


2.50 ± 0.11 ± 0.12 
2.56 ± 0.16 ± 0.13 
2.53 ± 0.09 ± 0.11 


0.45/3 (93.0) 
0.39/3 (94.2) 
0.95/7 (99.6) 


D°a; 


2.57 ± 0.11 ± 0.14 


3.19/3 (36.3) 


D077'(7r7rr?(77)) 
D°-n' (comb.) 


1.37 ± 0.14 ± 0.07 
1.73 ± 0.28 ± 0.08 
1.48 ± 0.13 ± 0.07 


0.05/3 (99.7) 
0.27/3 (96.6) 
1.55/7 (98.1) 


75*%° 


3.05 ± 0.14 ± 0.28 


4.73/7 (69.3) 


S*%(77) 
D*°r] (comb.) 


2.77 ± 0.16 ± 0.25 
2.40 ± 0.25 ± 0.21 
2.69 ± 0.14 ± 0.23 


4.20/7 (75.6) 

3.81/6 (70.2) 

10.48/14 (72.6) 


D*''u 


4.55 ± 0.24 ± 0.39 


4.05/7 (77.4) 


D* f)' {■WK'q^'y'y)) 
D'°r)' (comb.) 


1.37 ± 0.23 ± 0.13 
1.81 ± 0.42 ± 0.16 
1.48 ± 0.22 ± 0.13 


2.30/4 (68.1) 
0.68/2 (71.2) 
3.78/7 (80.5) 



V. RESULTS FOR THE B MEASUREMENTS 

The branching fractions measured in the different sec- 
ondary decay channels reconstructed in this analysis are 
given in Table IIIII (for missing entries in the Table; see 
the discussion on discarded sub- modes in Sec. IIII lT6|) . 
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FIG. 5: B{B'^ — >■ D^hP) (xlO"*) for the individual reconstructed D° and h^ decay channels (blue points) together with the 
BZ/ L^i? combination of this paper measurements (vertical yellow bands and the red points). The previous experimental results 
from BABAR 0, [lJ|, Belle |lll.[ia|. and CLEO [ly| are also shown (black points). The horizontal bars represent the statistical 
contribution alone and the quadratic sum of the statistical and systematic uncertainty contributions. The width of the vertical 
yellow band corresponds to ±1 cr of the combined measurement, where the statistical and systematic uncertainties are summed 
in quadrature. 
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FIG. 6: B{B° -^ D*°h'^) (xlO"'') for the individual reconstructed D° , D*°, and h° decay channels together with the BLUE 
combination of this paper measurements (vertical yellow bands and the red points). The blue squares (triangles) are for 
measurements with the sub-decay D*° — >■ D^tt" {D'^^j). The previous experimental results from BABAR 0, [lJ|, Belle |lll.[l2|. 
and CLEO [10| are also shown (black points). The horizontal bars represent the statistical contribution alone and the quadratic 
sum of the statistical and systematic uncertainty contributions. The width of the vertical yellow band corresponds to ±1 cr of 
the combined measurement, where the statistical and systematic uncertainties are summed in quadrature. 
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These branching fractions are combined using the so- 
called Best Linear Unbiased Estimate (BLUE) tech- 
nique [35| , that accounts for the correlation between the 
various modes. In the BLUE method the average value 
is a linear combination of the individual measurements, 



B = ^(a, • B., 



(13) 



where each coefficient at is a constant weight, not nec- 
essarily positive, for a given measurement Bi. The con- 
dition X]i=i cti — ^ ensures that the method is unbiased. 
The set of coefficients a = (ai, a2, ■•., cxt) is calculated 
so that the variance of branching fraction is minimal. 



E-^U 
U^E-^U' 



(14) 



where f/ is a i-component vector whose elements are all 
1 {U'^ is its transpose) and E is the (t x t) covariance 
matrix. The variance of branching fraction is then given 
by 



a^ Ea. 



(15) 



The covariance matrix E is evaluated for each source of 
systematics. Its matrix elements arc, for two modes i and 



Ej. 



Pij<^i 



(16) 



where Ui and Uj are the systematic uncertainties for the 
modes i and j, and pij is their correlation coefficient. 
We distinguish several types of systematic uncertainties 
according to their correlations between the modes: 

• full correlation, \pij\ ~ 1: neutrals (but uncertain- 
ties for 7r° and single 7 are independent), PID, 
tracking, number of BE, B{D*°), D*^uj polariza- 
tion in that mode, 

• medium correlation: B{D^), B(hP), whose correla- 
tions are taken from the PDG [2g| and range from 
2% to 100%, D^*^°p~ background in B° -> D^*^°tt^, 

• negligible correlation, \pij\ ~ 0: statistical uncer- 
tainties, PDF systematics, selection of intermediate 
particles, MC statistics. 

The total covariance matrix E is then the sum of the co- 
variance matrices for each source of uncertainty, plus the 
covariance matrix associated to statistical uncertainties. 
The systematic (statistical) uncertainty on the combined 
value of branching fraction computed by using Eq. (|15p 
where the error matrix includes only the systematic (sta- 
tistical) uncertainties. 

The combined branching fractions in data are given 
in Table IIVI with the x^ of the combination, the num- 
ber of degrees of freedom of the combination (ndof ) , and 
the corresponding probability (p- value). The individual 
branching fractions together with the combined value arc 



displayed in Figs. [5] and |6] and they are compared to the 
previous measurements by CLEO [ly|, BaBAR [3,[1J], and 
Belle [H [11. 

The results of this analysis, based on a data sam- 
ple of 454x10^ BB pairs, are fully compatible with our 
previous measurements [7|, |lj|, and also with those of 
CLEO [i3]- They are compatible with the measure- 
ments by Belle [111, [i2| for most of the modes, except for 
B^ -^ D(*)°?7, D*°a;, and £)*°7r°, where our results are 
larger. Those four branching fractions are from 2.5 to 
3.7 standard deviations (including systematic uncertain- 
ties) away from the latest measurements by Belle |lll.[l3|. 
Our measurements are the most precise determinations 
of the B{B° -^ £)(*)°ft,°) from a single experiment. They 
represent significant improvements with respect to the 
accuracy of the existing PDG averages [2a] . 

As a cross check we also perform the branching frac- 
tions measurements with the sub-data set of 88x 10^ BB 
pairs that we previously studied [7| . When studying the 
same decay modes we find values compatible with that of 
Ref. [3] with both statistical and systematic uncertainties 
lowered by significant amounts. In addition to the bene- 
fit from improved procedures to reconstruct and analyze 
the data, this updated analysis incorporates new decay 
modes, higher signal efficiency, and better background 
rejection and modeling. We use additional control data 
samples and measure directly in the data the relative ra- 
tio of the B^ — ?► D'^*^'^ p^ backgrounds. This analysis 
employs better fitting techniques and uses more sophisti- 
cated methods to combine the results obtained with the 
various sub-decay modes. 



VI. POLARIZATION OF B° -^ D*°lu 

The polarization of the vector-vector (VV) decay 
B'^ — > D*^u) has never been measured. Until now, it was 
supposed to be similar to that of the decay B~ — >• D*^p~ 
based on HQET and factorization- based arguments [3q . 
The angular distributions for the decay JS" — >■ D*'^lu are 
described by three helicity amplitudes: the longitudinal 
Hq amplitude and the transverse iJ+ and H„ amplitudes. 
In the factorization description of i? ^- VV decays, the 
longitudinal component Hq is expected to be dominant, 
leading to the fraction of longitudinal polarization, de- 
fined as 



JL ^ — — 



\Ho\ 



\Ho\ 



\H+\^ + \H- 



(17) 



predicted to be close to unity [31 I37H39J . 

Significant transverse polarizations were measured in 
B — )- (/)K* (see the review in (2611 ') and investigated as 
possible signs of New Physics [401 > but could also be the 
result of non-factorizable QCD effects [4l|. Similar ef- 
fects were studied in the context of SCET |22l, and are 
expected to arise in the B^ — )■ D*'^uj decay, in particular 
through enhanced electromagnetic penguin decays [42| . 
leading to significative deviation of /^ from unity. It has 



also been argued in SCET studies that non-trivial long 
distance contributions to the B^ — !• D*^uj amplitude may 
allow a significant amount of transverse polarization of 
similar size to the longitudinal polarization, leading to a 
value Jl ^ 0.5. 

Apart from the motivation of these phenomenological 
questions, the uncertainty on the angular polarization 
of B° — !• D*^uj affects the kinematic acceptance of this 
decay channel and therefore would be the dominant con- 
tribution to the systematic effects for its B measurement. 
Hence we measure the fraction of longitudinal polariza- 
tion for this decay mode. The analysis is performed with 



J 
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i?° -^ D*^uj candidates selected with the same require- 
ments as for the B analysis described in the previous 
sections. We consider the sub-decays D*^ -^ D'^-k^ and 
L)° -^ A'"7r+, A'"7r+7r°, Jir~7r+7r~7r+, and K^-k+t:^ . 



A. Description of the method 



The differential decay rate of B^ 

0^0 



decay D*^ -^ D^'tt" is 



D*"uj for the sub- 



d cos Oo'd cost 



^ 



[\H+\^ + \H-\^ + 2{Re{H+H*_)cos2x- Im{H+H*_)sm2x)] sin'^ Od' sin^^^-h (1^) 

{Re{H+H^ + H_H^) cosx - Im{H+H^ - H^H^) sin x) sin26'c. sin26l„, 



r 



where 6d* {0^) is the helicity angle of the D* {uj) me- 
son (see Sec. IIII(j3] for definitions). The angle x, called 
the azimuthal angle, is the angle between the D*^ and 
w decay planes in the B^ frame. Since the acceptance 
is nearly independent of Xi one can integrate over x to 
obtain a simplified expression: 



d^r 



■ ex 



A\Ho\^cos^eD'Cos^eu, (19) 



dcosOo'dcosO^ 

H\H+\^ + |H_ n sin^ Od' sin2 9^. 

This differential decay width is proportional to 

4/l cos^ Od' cos^ 61a; + (1 - /l) sin^ Od- sin^ O^^, (20) 

which is the weighted sum of purely longitudinal (/^ = 1) 
and purely transverse (/^ = 0) contributions. 

We employ high statistics MC simulations of exclusive 
signal samples of B^ -^ D*^lli decays with the two ex- 
treme configurations /l = and 1 to estimate the ratio 
of signal acceptance, £o/eii of /l = events to /l = 1 
events. The longitudinal fraction /i, can be expressed in 
terms of the fraction of background events (7) and the 
fraction of /l = 1 events in the observed data sample 
(a): 



II 



a -t- (1 — a — 7) • — 
£1 



(21) 



The fraction 7 is taken from the fit of /S.E for a sig- 
nal region |A£;| < 2.5 ctab and mEs > 5.27 GeV/c^, 
where a^E is the fitted A£' width of the signal distri- 
bution, ranging from 20.8 to 23.3 MeV depending on the 
mode. The fraction a is determined from a simultane- 
ous 2-dimensional fit to the distributions of the helicity 



angles cos 9^ and cos Oij, , for B'^ -^ D*°uj candidates se- 
lected in the same signal region. The correlation between 
cos9i^ and cos9d' is found to be negligible. 



The signal shapes are described with parabolas (see 
Eq. (PU]) ). except for the cos9i^ distribution of /^ = 
signal events, which is described by a non-parametric 
PDF based on the MC simulation. As the signal dis- 
tribution of cos9i^ is distorted around zero because of 
the selection cut on pion momentum and on the co boost 
(see Sec. IIIIB3|) . The signal PDF parameters are fixed 
to those fitted on the D*'^lu simulations. The shape 
of the cos 9^ and cos 9d* background distributions is 
taken from the data sideband \AE\ < 280 MeV and 
5.235 < mEs < 5.270 GeV/c^. The consistency of the 
background shape was checked and validated for various 
regions of the sidebands in data and generic MC simula- 
tions. Possible biases on /^ from the fit are investigated 
with pseudo-experiment studies for various values of /l 
from 5 to 95 %. No significant biases are observed. 



An additional study is performed with an embedded 
signal MC simulation, i.e. with signal events modeled 
from various different fully simulated signal samples and 
with a generated value /^ ~ 90% (as expected from 
HQET [si, El]). A small bias on the fitted value of 
/l is observed (~ 14% of the statistical uncertainty). 
This bias is due to the slight difference on the descrip- 
tion of the signal shape for the cos 9ui distribution and 
for /i = 0, modeled by a non-parametric PDF, and that 
of the actual shape obtained from the embedded signal 
MC simulation. This bias is corrected later on and we 
assign a systematic uncertainty on that correction. 
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FIG. 7: Fitted distributions of the helicity cos So* and cos 6^ in the decay channel B" — >■ D*'^lu for the D*^ decay modes K~-k^ 
(a) and (b), K~-k'^tt^'' (c) and (d), K~n^Tr~-jv^ (e) and (f),and KsTv'^tv~ (g) and (f). The dots with error bars are data, 
the curves are the various PDF contributions: the sohd blue (blue) is the total PDF, the dash-dot (grey) is the background 
contribution, the long dash (blue) is the /l = 1 signal part and the dots (red) is the /l = signal. 



B. Statistical and systematic uncertainties 

The statistical uncertainty on /^ is estimated with 
a conservative approach by varying independently the 
values of the two fitted parameters a and 7 by ±1 a 
in Eq. (pij) . An extended study based on MC pseudo- 
experiments accounting for the correlations between a 
and 7 gave slightly smaller uncertainty. 

The uncertainty due to the signal shape in the simul- 
taneous 2-dimcnsional fit to cos 9^ and cos Oo* is mea- 
sured using the control sample B^ -^ £)*°7r~, with 
£)*o ^ £,0^0 g^j^j £)0 ^ K-7r+. This mode was cho- 
sen for its high purity and for its longitudinal fraction 
/l = 1, which enables us to directly compare its shape 
to our signal /l = 1 . The distribution of the helicity an- 
gle of the D*'^ is found to be wider in the data than in the 
MC, this difFerencc being parameterized by a parabola. 
The uncertainty on the signal shape is then measured by 
refitting a, with the signal PDF being multiplied by the 
correction parabola. The relative difference is then taken 
as the uncertainty. 

We assign a systematic effect due to the correction that 



we apply for the observed small bias on the value for /l, 
when fitting the embedded signal MC simulation (see the 
discussion at the end of Sect. IVI Ap . 

The uncertainty due to the background shape is mea- 
sured by refitting a with the background shape fitted in 
a lower data sideband |A£:| < 280 MeV and 5.200 < 
Toes < 5.235 GeV/c^. The relative difference is then 
taken as the uncertainty. 

An uncertainty is assigned to /^ due to the assumption 
of the acceptance being independent of x- The accep- 
tance of the MC simulation signal is measured in bins of 
X and fitted with a Fourier series to account for any devi- 
ation from fiatness. The resulting fitted function is used 
as a parametrization of the acceptance dependency to x 
and multiphed to the decay rate (as written in Eq. (fT9)) ). 
We then perform in a study based on pseudo-MC simu- 
lation experiments where the events are generated from 
this new decay rate. The resulting cos 9^ x cosBd* dis- 
tributions are fitted with the procedure described above. 
A small bias is observed and its value is assigned as a 
systematic effect. 

The uncertainty on the efficiency ratio £o/ei, from the 
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TABLE V: Total relative uncertainties computed in data on 
the measurement of /^ in the decay channel B" — >■ D'^'lj, 
with D*° -> D°7v° and D° -^ K-n+, K'Tr+n", K--k+-k--k+, 
and K%-k'^-k~ . 



TABLE VL Values of a fitted in data, of the background 
fraction 7 and of the acceptance ratio eo/ei, with the corre- 
sponding values of the longitudinal fraction /l after the bias 
correction. The first quoted uncertainty is statistical and the 
second systematic. 



Sources 



Ml/Il (%) 
K-K K3tt Ktvk° K'-Itttv 



Signal PDFs 
Correction of the bias 
Background PDF 
Limited MC statistics 
Flat acceptance vs. x 



Total systematic 


3.1 


5.6 


4.8 


8.6 


Statistical uncertainty 


9.6 


16.3 


16.3 


25.6 


Total uncertainty 


10.0 


17.2 


17.0 


27.0 



limited amount of MC statistics available, is calculated 
assuming Eq and £1 to be uncorrclated, while individual 
uncertainties on eg and ei are calculated assuming a bi- 
nomial distribution. 

The various relative uncertainties are displayed in Ta- 
ble |V] for the data and are found to be compatible with 
those calculated in MC simulations. The dominant un- 
certainty is statistical. Among the various systematic 
sources, the largest contribution comes from the signal 
and background parametrizations, and for some modes 
on the assumption of a flat acceptance versus x- 

As a check, the Jl measurement is applied in data 
first on the high purity and high statistics control sam- 
ple B- ~> D 



*o. 



with D*" -> D"n" and D'' -> K-jt 



This decay channel is longitudinally polarized, i.e. Jl = 
1. The fit of cos 9d* in data yields a value of /l compat- 
ible with one, as expected. 



D° mode 



(%) 7 (%) £o/£i 



h {%) 



2.5 


2.9 


2.4 
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33.4±2.7 52.0±1.9 
18.8±2.3 71.2±2.5 
19.6±2.1 76.0±2.3 
24.9±4.2 66.0±4.9 



1.093±0.012 64.8± 6.5±2.1 
1.068±0.017 60.8±10.3±3.6 
1.109±0.021 76.9±13.0±3.8 
1.092±0.016 66.7±18.3±6.2 



Combi. 



fL = (66.5 ±4.7 ±1.5)% 




D°(^ K'k^k"] 
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D°(^ K"it*) 
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f L (%) for b" -^ D'°m 

FIG. 8: Measurements of /z, with the four D^ modes in data. 
The yellow band represents the BL UE combination. 



C. Results for the fraction of longitudinal 
polarization Jl 

The fitted data distributions of the cosine of the helic- 
ity angles are given in Fig. [T] The measurements for each 
D^ decay channel are then combined with the BL UE sta- 
tistical method [H] (see Sec. |V| with xV^dof = 1.01/3 
[i.e.: a probability of 79.9%). The measured values of 
Jl, ol, 7 and £o/£i are given with the details of the com- 
bination in Table I VII and in Fig. [H The final result is 
Jl = (66.5 ± 4.7 ± 1.5)%, where the first uncertainty is 
statistical and the second systematics. This is the first 
measurement of the longitudinal fraction of B^ — > D*'^oj, 
with a relative precision of 7.4%. 

This value differs significantly from the HQET pre- 
diction Jl = (89.5 ± 1.9)% Jsi, 13. This significant 
transverse amplitude in the B^ — > D*^uj decay channel 



may arise from the same mechanism as the one that is 
responsible for the transverse polarization observed in 
B -^ <f>K* . It however supports the existence of effects 
from non-trivial long distance contributions to the decay 
amplitude of S° -^ D*°lo as predicted by SCET stud- 
ies l2a. 



VII. DISCUSSION 

A. Isospin analysis 

The isospin symmetry relates the amplitudes of the de- 
cays B- -^ D(*)07r-, fiO -^ D'^*^+n- and 5° -^ D'^*^°tt°, 
which can be written as linear combinations of the isospin 
eigenstates Aj ^o, I = 1/2, 3/2 0,1451, 
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y^(Z?(*)"^-) 


— V 3.43/2, _D(*) : 




(22) 


A{D'^*^+'K-) 


= l/V3.43/2,£i( 


) + V 2/3.4i/2,D(* 


): 


^p(*)o^O) , 


= v/2/3-43/2,_D( 


.) - VV3-4i/2,D( 


)i 


eading to: 








yl(I?(*)"7r-) = 


A{D^*^ + TT-) + 


^/2y^p(*^V"). 


(23) 



The relative strong phase between the amphtudes 
-^1/2, D(*) ^-nd -^3/2, _D(*) is denoted as 5 for the Dtt sys- 
tem and 5* for the D*7r system. Final state interactions 
between the states Z)^*-'°7r° and D^*'^tt^ may lead to a 
value of J^*-* different from zero and, through construc- 
tive interference, to a larger value of B for £)(*'°7r° than 
the prediction obtained within the factorization approx- 
imation. One can also define the amplitude ratio R^*\ 
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1-^1/2, D(') I 
V2|.43/2,D(.) I 



In the heavy-quark limit, the factorization model pre- 
dicts 0, 113 5^*^ = 0{AQCD/mb) and i?(*) = 1 4- 
©(ylqCD/'TT-b), where nib represents the b quark mass and 
where the correction to "1" is also suppressed by a power 
of l/Nc, with TVc the number of colors. SCET pfl-ig pre- 
dicts that the strong phases 5^*'> (ratios R^*') have the 
same value in the D-k and I?*7r systems and significantly 
differ from (1). 

The strong phase 5^*' can be computed with an isospin 
analysis of the D'^*'^it system. We use the world av- 
erage values provided by the PDG [2g for B{B~ — > 
D(*)°7r-), B{B° -^ £)(*)+7r-) and for the B lifetime ratio 
t{B+)/t{B°). The values oiB{B° -^ D'^*^°n°) are taken 
from this analysis. We calculate the values of iJ^*-* and 
R^*^ using a frequcntist approach |48j : 



(24) FIG. 9: Combined ratios B{B° -^ D'-°h^)/B{B^ -^ D"h°) 
measured in tills paper compared to theoretical prediction by 
SCET [221 (vertical solid line). The vertical band represent 
the estimated theoretical uncertainty from SCET (for the case 
where h" is the uj meson see text). 



6 = (29.0+ti) 
for Dtt final states, and 



R - {69.2tli) %, 



i29.5tllr 



R* = (67.0 



-l-4.8\ 
4.7) 



%, 



(25) 



(26) 



for D*Tr final states. 

In both Dtt and D*tt cases, the amplitude ratio is 
significantly different froin the factorization prediction 
jj(*) _ I -pjjg strong phases are also significantly dif- 
ferent from zero and are equal in the two systems Dtt 
and D*n (0° is excluded at 99.998% and 99.750% of con- 
fidence level, respectively), which points out that non- 
factorizablc FSI arc indeed not negligible. Those results 
confirm the SCET predictions. 



B. Comparison to theoretical predictions for 

Table Eni compares the B{B" -^ D^*^"h°) measured 
with this analysis to the predictions by factorization [3|, 



II 



1*0 



0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 
BF(B° -^ X ri')/ BF(B° ^ X ri) 



FIG. 10j_ Combined ratios B{B° ^ D*'''rj')/B(B° -^ D''°-q) 
and B{B° -i> D°-q')/B{B^^ -> D^-q) measured in this paper 
compared to theoretical prediction by SCET [22| (vertical 
line) and from factorization [4y| (vertical bands). 



Hliili^] and pQCD [Tl,[l|. We confirm the conclusion 
by the previous BaBAR analysis [3] : the values measured 
are higher by a factor of about three to five than the 
values predicted by factorization. The pQCD predictions 
are closer to experimental values but are globally higher, 
except for the _D(*)''7r" modes. 
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TABLE VII: Comparison of the measured branching fractions B, with the predictions by factorization [J, [la . |49| . |50| and 
pQCD [la.[l9l|. The first quoted uncertainty is statistical and the second is systematic. 



B{B° ->■) (xlO"'') This measurement 



Factorization 



pQCD 



,0^0 


2.69 ± 0.09 ± 0.13 


0.58 yj; 0.70 [3] 




2.3-2.6 


,.0^0 


3.05 ± 0.14 ± 0.28 


0.65 [16J; 1.00 [3] 




2.7-2.9 


,0^ 


2.53 ± 0.09 ± 0.11 


0.34 [16]; 0.50 3] 




2.4-3.2 


„0^ 


2.69 ± 0.14 ± 0.23 


0.60 [3] 




2.8-3.8 


<0i^ 


2.57 ± 0.11 ± 0.14 


0.66 [16]; 0.70 [3] 




5.0-5.6 


"0^ 


4.55 ± 0.24 ± 0.39 


1.70 [3] 




4.9-5.8 


firi' 


1.48 ± 0.13 ± 0.07 


0.30-0.32 [50|; 1.70-3.30 


[49] 


1.7-2.6 


,.0^^, 


1.48 ± 0.22 ± 0.13 


0.41-0.47 [49.1 




2.0-3.2 



The ratios of the B are given in Table I Villi It should 
be noted that the values of these ratios are not computed 
directly from those quoted in Table llVi as we take ad- 
vantage of the fact that common systematic uncertainties 
cancel between ZJ^/i^ and D*^h^ modes. Therefore the 
ratios of the B are first calculated for each sub-decays of 
D^ and /i°, and then combined with the BLUE method. 
The ratios B{B^ -> D*°h^)/B{B° ^ D^h°) for h^ = tt", 
?7, and if are compatible with 1. All are displayed in 
Fig. [9] together with the theoretical predictions. 

Factorization predicts the ratio B{B^ — >■ 
D'^*)Ori')/B{B° -^ D^*^°r]) to have a value between 
0.64 and 0.68 [i^l, related to the 77 — 77' mixing. 
Those ratios are also given in Table IVIIII and Fig. [TU] 
compares the theoretical predictions with our ex- 
perimental measurements. The measured ratios are 
smaller than the predictions and arc compatible 
at the level of less than two standard deviations. 
The SCET gives also a prediction about the ratio 
B{B^ ^ D^*^"r]')/BiB° -^ L>(*)o?7) ~ 0.67, which is 
similar to the prediction by factorization. 

SCET f20l - [2^ does not predict the absolute_value 
of the B but it predicts that the ratios B{B'^ — > 
D*^h°)/B(B^ -> D°/i°) are about equal to one for h^ = 
TT*^, 77 and 77'. For h^ — u that prediction holds only for 
the longitudinal component of i?° — s- D*'^ll!, as non trivial 
long-distance QCD interactions may increase the trans- 
verse amplitude. We measure the fraction of longitudinal 
polarization to be Jl = (66.5 ± 4.7(stat.) ± 1.5(syst.))% 
in the decay mode i?° — > D*'^uj, and find that the ra- 
tio B{B° -> D*°uj)/B{B" -^ D°uj)Js significantly higher 
than one, as expected by SCET [23 |. 



VIII. CONCLUSIONS 

We measure the branching fractions of the color- 
suppressed decays B° — > D'-*^^h°, where h° = 7r°, 77, 
uj, and 77' with 454xl0^i3i3 pairs. The measurements 
are mostly in agreement with the previous results [3, Hol - 
[l2l . [Til and are the most precise determinations of the 
B{B^ — )■ D^*^^h^) from a single experiment. They repre- 



TABLE VIII: Ratios of branching fractions S(B<' -^ D*°h°)/ 
Z3(S° ^ D"h°) and B{B'> -, D^*^'>ri')/B{B° ^ £>(*)"r;). The 
first uncertainty is statistical, the second is systematic. 



B ratio 



This measurement 



D*0tt"/D07t0 
D*0r,(77)/D0,7(77) 
D*°ri(TnnT°)/D"r){TnT7T^'') 
D*°r]/D"r) (Combined) 

D*°LU/D°LU 

D*°ri' {■KTxr])/D°r]' {-K-wq) 
D'^^-q' /D°'q' (Combined) 



1.14 ± 0.07 ± 0.08 
1.09 ± 0.09 ± 0.08 
0.87 ± 0.12 ± 0.05 
1.03 ± 0.07 ± 0.07 
1.80 ± 0.13 ± 0.13 

1.03 ± 0.22 ± 0.07 
1.06 ± 0.38 ± 0.09 

1.04 ± 0.19 ± 0.07 



D*"ri'/D*Ori 



0.54 ± 0.07 ± 0.01 
0.61 ± 0.14 ± 0.02 



sent significant improvements with respect to the accu- 
racy of the existing PDG averages |26l |. 

For the first time we also measure the fraction of longi- 
tudinal polarization /^ in the decay mode B° — >■ D*'^u> to 
be significantly smaller than 1. This reinforces the con- 
clusion drawn from the B measurements on the validity 
of factorisation in color-suppressed decays and supports 
expectations from SCET. 

We confirm the significant differences from theoretical 
predictions by factorization and provide strong con- 
straints on the models of color-suppressed decays. In 
particular our results support most of the predictions of 
SCET on fiO -^ D^*^^h° ^^^. 
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• Ap is the value at the maximum of the function, 

• Xp is the peak position, 

• ap is the width of the peak defined as the width at 
half-height divided by 2^2 In 2 ~ 2.35, 

• r is an asymmetry parameter. 

• Pi/2 described the width of the respective tails. 



=F 1 . Left 



The positions 3; 1,2 are Xp + (Tp\/2\n2 
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